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Abstract. We derive photometric, structural and dynamical evolution-related parameters of five as yet unstudied 
low-contrast open clusters located in the third quadrant using 2MASS data. The target clusters are Czernik31, 
Czernik32, Haflher9, Haffher 11 and Trumpler 13. We apply a statistical field-star decontamination procedure to 
infer on the intrinsic colour-magnitude diagram (CMD) morphology which is critical for such low-contrast objects. 
Consequently, it became possible to derive accurate reddening, age, distance from the Sun and Galactocentric 
distance for the five clusters. In the structural and luminosity /mass-function analyses we apply a colour- magnitude 
filter which encompasses the cluster evolutionary CMD sequences and excludes stars with discrepant colours. Using 
this procedure we derive core and limiting radii, mass function (MF) slope, total mass, mass density and relaxation 
time. We derive ages in the range 140 - 1 100 Myr, Galactocentric distances within 7.7 - 11.4 kpc, and total masses 
within 360-2 900 Mq. Reflecting large-scale mass segregation, the MF slope in the core is significantly flatter than 
that in the halo of the five clusters. Although some of the present clusters are relatively younger than the Gyr-old 
clusters, they present evidence of advanced dynamical evolution. This kind of study has become possible because 
of the photometric uniformity and spatial coverage of 2MASS which allows a proper subtraction of the field-star 
contamination on the target CMDs. The present study indicates that low-contrast clusters can be studied with 
2MASS, particularly after field-star subtraction, which is important since most of the unstudied open clusters 
belong to this class. 

Key words. (Galaxy:) open clusters and associations: general; Galaxy: structure 
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covery and analysis of two faint open clusters in Cygnus 
(Bica, Bonatto, & Dutra 2004) and three in other parts 
of the sky (Bonatto, Bica & Dutra I2UM|> using 2MASS 
photometry. In Bonatto & Bica (2005) we used 2MASS 
to obtain a homogeneous set of parameters, compare dy- 
namical states and derive mass functions (MFs) for 11 
relatively populous, nearby open clusters. 

Considering the above we decided to study five open 
clusters in the third quadrant using 2MASS data and the 
techniques outlined in Bonatto & Bica (2005). Because 
the present sample is composed of low-contrast clusters 
(Sects. EJand^J) we include in the analysis an algorithm to 
statistically decontaminate the observed CMDs of Galaxy 
(field) stars to better define the CMD morphology. This 
step is essential to unambiguously separate low-contrast 
open clusters from fluctuations in field-star counts. In ad- 
dition, colour-magnitude filters encompassing the cluster 
evolutionary CMD sequences are subsequently used in the 
structural and luminosity/mass-function analyses. 

For the present clusters we derive the age, reddening, 
distance from the Sun, Galactocentric distance, core and 
limiting radii, mass and density, MF slope and relaxation 
time. We derive parameters for the cluster as a whole as 
well as the core and halo subsystems. Besides the interest 
in obtaining parameters of as yet unstudied clusters, this 
paper represents a test of the limits of 2MASS (and our 
techniques) when dealing with more distant, low-contrast, 
low-Galactic latitude clusters. 

Uncertainties throughout this paper correspond to ler 
Poisson statistics. 

This paper is organized as follows. In Sect.|21we present 
the sample and show XDSS R images of the clusters. In 
Sect. 03 we present the 2MASS photometric parameters 
and introduce the field-star decontamination procedure. 
In Sect.^lwe analyze the spatial structure of the clusters. 
In Sect. El we derive the mass functions and discuss stel- 
lar content properties. In Sect. we compare the present 
clusters with nearby, populous open clusters which span a 
range of dynamical states. Concluding remarks are given 
in Sect.d 



2. The third quadrant (3Q) open cluster sample 

The clusters selected for the present study are Czernik31 
(OC1-625, ESO560SC3), Czernik32 (Ki24, OC1-683, 
BH11, ESO494SC20), Haffner9 (OC1-600), Haffnerll 
(OC1-657, BH3, ES0429SC3, and Trumplerl3 (Cr219, 
OC1-815, BH94, ES0127SC17). Additional desig nation s 
above are from Alter, Ruprecht & Vanisek l|1970|l . 
Lauberts l|1982|) . and van den Bergh & Hagen (|1975|l . In 
what follows we will use the acronyms Cz, Haf and Tr to 
refer to the clusters. For simplicity we will refer to the 
present clusters as the 3Q sample. 

In Figs. ^ and |2] we present XDSS R images of the 
clusters in the 3Q sample. DSS and XDSS images can 
be extracted from the Canadian Astronomy Data Centre 



(CADC 3 ). In all panels a central concentration of stars 
with varying star-density contrast with respect to the 
background can be seen. 

Accurate central coordinates are fundamental for the 
photometric, and primarily structural, analysis of open 
clusters. This is particularly critical in the derivation of 
the core radius and mass of small clusters (see e.g. Bonatto 
& Bica 2005 ). For the 3Q sample clusters we indicate in 
Tabler^diffcrent sources of coordinates. The original Alter, 
Ruprecht & Vanisek (1970) coordinates (1st line for each 
cluster in Tabled are not accurate enough for the present 
purposes of deriving core, halo and overall parameters. 
Consequently we measured in the DSS and XDSS images 
more precise optical centers for each cluster. Figs. H and 
El are centered on these coordinates. 

3. The 2MASS photometry and near-infrared 
CMDs 

The VizieR 4 tool was used to obtain J, H and K s 2MASS 
photometry. The central coordinates corresponding to the 
density peaks in the 2MASS photometry of each cluster 
were derived by examining histograms of the number of 
stars in 0.5' bins of right ascension and declination. In 
what follows we refer to the coordinates that maximize 
the density of stars (Sect.0J at r = 0' as the cluster center. 
These central coordinates are given in columns 2 and 3 of 
Tabled In a few cases the central coordinates are slightly 
shifted with respect to the optical centers (Figs. 2] and 
F2J). For each cluster we extracted photometry of the stars 
contained in a circular area with radius R ex tr (column 7 
of TableQ centered on the respective central coordinates. 

Because of the relatively low Galactic latitudes of the 
3Q sample clusters (column 5 of Tabled, the target fields 
present significant field-star contamination, mostly from 
disk stars. The relative amount of this contamination with 
respect to the central field turns out to be important in 
the identification and analysis of small and/or poor open 
clusters. To illustrate the relative contribution of the field 
stars with respect to varying spatial areas we show in the 
left panels of Figs. to the J x (J - H) CMDs of the tar- 
get clusters in two different extractions. The correspond- 
ing (same area) offset fields are shown in the middle pan- 
els. The smaller extraction (panels (a)) is that maximizing 
the CMD density contrast with respect to the offset field. 
The smaller radius corresponds to about 4x the core ra- 
dius (Rcoro, Sect. HJ. The radius of the larger extraction 
(panels (b)) is intermediate between that of the smaller 
extraction and the limiting radius (Rn m , Sect.0}. The off- 
set fields were built with the stars in rings more external 
than the limiting radius. 

The cluster CMDs contrast in morphology and den- 
sity with respect to the offset fields, although field-star 
contamination in some cases is important, e.g. Cz31 and 
Cz32. To minimize the visual effect of the field stars on 

3 http://cadcwww.dao.nrc.ca/ 

4 \http:/ /vizier. u-strasbg.fr/viz-bin/VizieR ?-source=II/24.6\ 
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Fig. 1. XDSS R images. Left panel: Haf9 (6' x 6'). Middle panel: Haf 11 (7' x T). Right panel: Cz32 (6' x 6' 




Fig. 2. XDSS R images. Left panel: Cz31 (7' x T). Right panel: Tr 13 (6' x 6'). 



the cluster CMDs we applied a decontamination proce- 
dure which considers two different approaches for field-star 
subtraction: (i) absolute counts and (ii) nearest-neighbour 
density counts. In case (i) the algorithm divides the CMD 
in boxes of dimensions AJ = 0.75 and A(J — H) = 0.075, 
both for the cluster and offset field. The effective number 
of cluster stars in a given CMD box corresponds to the ab- 
solute difference in the number of cluster and offset-field 
stars. After this step the algorithm randomly excludes 
from each cluster CMD box the required number of stars; 
boxes resulting with a negative number of stars are left 
blank. In case (ii) we calculate the star density in CMD 
boxes with dimensions AJ = 0.35 and A(J — H) = 0.035. 
The observed star density in a given box corresponds to 
the average number of stars in it and the 1st, 2nd and 3rd- 
neighbouring row and columns, both for the cluster and 
offset-field regions. Subtraction of the observed offset-field 
density from that of the cluster yields the effective cluster 
star density in each CMD box. Multiplying the ratio of 
the effective/observed star density by the observed num- 
ber of stars yields the effective number of cluster stars in 
a given CMD box. From this point on the algorithm pro- 
ceeds as in approach (i) to exclude stars in each CMD 
box. Comparison of the resulting decontaminated CMDs 



with a variety of offset fields indicated that the best re- 
sults were obtained with approach (ii), probably because 
it is less sensitive to field-star count fluctuations than ap- 
proach (i). To maximize the statistical representativity of 
the field-star counts we consider as offset field the region 
1.5 x Ru m < r < R ex t- The ratio between the cluster and 
offset-field spatial areas is considered when we calculate 
the respective number (or density) of stars. Approach (i) 
is similar to that used by Kerber & Santiago (2005), and 
(ii) is similar to Mighell et al. pHHfij) . 

The resulting field-star decontaminated CMDs are 
shown in the right panels of Figs. to [7| The decon- 
taminated CMDs are cleaner than the observed ones (left 
panels) and present morphologies typical of open clusters 
of different ages. In particular the giant clumps of Haf 11 
(Fig.^J and Cz 32 (Fig.0 show up in the respective CMDs 
which in turn helps constrain their ages and distances from 
the Sun. The field-star decontaminated CMDs are used to 
better define the intrinsic CMD morphology and isochrone 
fit. 

Photometric parameters were derived by means of so- 
lar metallicity Padova isochrones (Girardi et al. |2002) 
computed with the 2MASS J, H and K s filters 5 , fitted to 

5 http:/ /pleiadi.pd. astro, it/ isoc_photsys. 01/isoc_photsys. 01 . htmi 
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Table 1. 3Q cluster coordinates and 2MASS photometric parameters 



Object 


a(J2000) 


<5(J2000) 


i 


b 


z 


R-cxtr 


S' c 


Ago 


E(B - V) 


(m - M)o 


Mj 


d© 


Scale 


dec 




(h:m:s) 


(":':") 


n 


(°) 


(PC) 


(') 




(Myr) 


(mag) 


(mag) 


(mag) 


(kpc) 


(PC/') 


(kpc) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


Haf 9+ 


07:24:30 


-17:00 


























Haf9 


07:24:42.0 


-17:00:10.0 


231.80 


-0.59 


-20 ± 2 


30 


3.9 ± 0.9 


140 ± 20 


0.54 ± 0.03 


11.4 ± 0.1 


3.1 


1.9 ± 0.2 


1.78 


9.3 ± 0.2 


Haf 11 + 


07:35:24 


-27:44 


























Hafll 


07:35:30.0 


-27:41:58.6 


242.39 


-3.52 


-322 ± 15 


40 


2.8 ± 0.5 


890 ± 150 


0.36 ± 0.03 


13.6 ± 0.1 


1.1 


5.2 ± 0.2 


0.66 


11.4 ± 0.2 


Cz 32^ 


07:50:24 


-29:50 


























Cz 32 


07:50:30.0 


-29:50:45.0 


245.89 


-1.73 


-121 ± 8 


30 


4.7 ± 1.2 


1120 ± 200 


0.66 ± 0.03 


13.0 ± 0.1 


1.4 


4.0 ± 0.2 


0.85 


10.3 ± 0.3 


Cz31 f 


07:36:48 


-20:30 


























Cz31 


07:36:59.0 


-20:30:35.0 


236.27 


+0.27 


+ 11 ± 1 


50 


3.4 ± 1.2 


180 ± 20 


0.06 ± 0.03 


11.7 ± 0.1 


3.2 


2.2 ± 0.2 


1.53 


9.4 + 0.2 


Tr 13 + 


10:23:48 


-60:05 


























Tr 13 


10:23:49.2 


-60:07:32.5 


285.51 


-2.34 


-100 ± 6 


15 


1.9 ± 0.4 


320 ± 40 


0.30 ± 0.03 


11.9 ± 0.1 


2.8 


2.4 ± 0.2 


1.41 


7.7 + 0.2 



Tabic Notes. Col. 6: distance to the Galactic plane. Col. 7: Extraction radius of the 2MASS photometry. Col. 8: Observed contrast parameter <5 C . 
Col. 9: Best-fit isochronc age. Col. 10: Colour-excess resulting from the isochrone fit. Col. 11: intrinsic distance modulus. Col. 12: faint-absolute 
magnitude limit of the stars in the radial-density profiles. Col. 13: distance from the Sun. Col. 14: parsec to arcmin scale. Col. 15: Galactocentric 
distance, (f): coordinates from Alter et al. 119701 . 



the field-star decontaminated CMDs of the central cluster 
regions in Figs. to □ (right panels). The 2MASS trans- 
mission filters produced isochrones very similar to the 
Johnson ones, with differences of at most 0.01 in (J — H) 
(Bonatto, Bica & Girardi 2004). For reddening and ab- 
sorption transformations we use Ry = 3.2, and the rela- 
tions Aj = 0.276 x A v and E(J - H) = 0.3 3 x E(B - V), 
according to Dutra, Santiago & Bica l|2002l and references 
therein). 2MASS photometric uncertainties as a function 
of magnitude are discussed in Soares & Bica t2002J and 
Bonatto, Bica & Santos Jr. <|2005[) . 

Taking into account the uncertainties associated with 
the isochrone fit we derive for each cluster in the 3Q sam- 
ple the age (column 9 in Table ^) colour excess E(J — H) 
and the corresponding E(B — V) (column 10), the absolute 
distance modulus (m — M)o (column 11) and the distance 
from the Sun d Q (column 13). The Galactocentric distance 
(column 15) was calculated using 8.0 kpc as the distance 
of the Sun to the center of the Galaxy (Reid I1993fl . The 
resulting isochrone fits are shown in the field-star decon- 
taminated CMDs (right panels of Figs. El to 0. 

3Q cluster ages are 140 ± 20 Myr (Haf 9), 180 ± 20 Myr 
(Cz31), 320 ±40 Myr (Trl3), 890 ±150 Myr (Hafll) and 
1.1 ± 0.2 Gyr (Cz31). Except for Cz31, the remaining 
clusters in the 3Q sample are considerably reddened with 
0.30 < E(B - V) < 0.66. They are distant from the Sun 
with 1.9 < d Q (kpc) < 5.2 and, except for Tr 13, the re- 
maining 3Q clusters are located outside the Solar circle. 

The relatively low-contrast nature of the clusters in 
the 3Q sample is apparent in the XDSS R images (Figs.^ 
and|2J) and is reflected in the low values of the observed 
contrast parameter S' c — 1.9 — 4.7 (column 8 of TableQ}. 
S' c represents the ratio of the number of stars in the clus- 
ter region with respect to that (same area) in the offset 
field. Additional clues to their nature as open clusters are 
(i) the cluster-like field-star decontaminated CMD mor- 
phology and isochrone fits (Figs. EltoEJ , (ii) the King-like 




Fig. 3. J x (J - H) CMDs of the central 4' extraction 
of Haf 9 (panel (a)) and respective same- area offset field 
(panel (b)). The corresponding field-star decontaminated 
CMD is in panel (c). Bottom panels: same as above for 
the 6' extraction. Solid line: best-fit Padova isochrone 
with age = 140 Myr. Dashed line: colour-magnitude filter 
(Sect.^J) used in the radial and luminosity/mass- function 
analyses. 
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Fig. 4. Same as Fig. for Hafll. The Padova isochrone Fig. 5. Same as Fig. [3] for Cz32. The Padova isochrone 
corresponds to the age 890 Myr. Stars above the turnoff corresponds to the age 1.1 Gyr. Stars above the turnoff 



may be blue-stragglers. 



may be blue-stragglers. 



radial-density profiles (Sect.^J, and (Hi) the resulting MFs 
(Sect. El). 

4. Cluster structure 

Structural parameters of the 3Q clusters were derived by 
means of the radial density profile, i.e. the projected num- 
ber of stars per area around the cluster center. Before 
counting stars we applied the colour-magnitude filters 
(dashed line in the right panels of Figs. 03 to |7J) to the 
observed CMDs to discard stars with discrepant colours 
(mostly Galactic field stars) and maximize cluster mem- 
bership. This filtering procedure was previously applied 
in the analysis of the open clusters M 67 (Bonatto & Bica 
HUD), NGC 188 (Bonatto, Bica & Santos Jr. 0005)1 and 
NGC3680 (Bonatto, Bica & Pavani 13004) . This proce- 
dure minimizes the probability of field-star contamination, 
spurious detections and the increase of photometric un- 
certainties at faint magnitudes. We give in column 12 of 
Tabled the resulting absolute faint- magnitude (Mj) limit 
of the stars included in each radial-density profile. 

The radial density profiles were obtained by counting 
stars inside concentric annuli with steps of 1.0' (Cz 31 and 
Hafll) and 0.5' (Cz 32, Haf 9 and Tr 13) in radius. In each 



case the field-star contribution level corresponds to the 
average number of stars included in the outermost ring. 

The resulting radial density profiles are shown in 
Fig. |5] For absolute comparison between clusters we scale 
the radius in the abscissa in parsecs, and the number den- 
sity of stars in the ordinate in stars pc -2 using the dis- 
tances derived in Sect. 01 The statistical significance of 
each profile is reflected in the ltr Poisson error bars. The 
cluster limiting radius (Rii m ) can be estimated considering 
the fluctuations in the radial density profile with respect to 
the field stars. In this sense, Rn m corresponds to the region 
where the density profile merges with the background. For 
regions beyond Rn m the null-contrast between cluster and 
field-star density would produce prohibitive Poisson errors 
and meaningless results. For practical purposes the bulk 
of the cluster stars are contained within Rn m . The tidal 
radius derived from the three-parameter King 1|1962|1 pro- 
file can only be obtained for much more populated open 
clusters, such as e.g. N GC 188, M 67, NGC 2477 and M 26 
(Bonatto & Bica ESfM)) . 

Cluster structural parameters were derived by fitting 
the two-parameter King fl966a! surface density profile 
to the background-subtracted radial distribution of stars. 
The two-parameter King model essentially describes the 
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Fig. 6. Same as Fig. [21 for Cz31. The Padova isochrone 
corresponds to the age 180 Myr. 
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(J-H) (J-H) (J-H) 

Fig. 7. Same as Fig. [21 for Tr 13. The Padova isochrone 
corresponds to the age 320 Myr. 



intermediate and central regions of normal clusters (King 
I1966bl Trager, King & Djorgovski 1995). The fit was per- 
formed using a nonlinear least-squares fit routine that uses 
the errors as weights. The best-fit solutions are shown in 
Fig. El superimposed on the respective radial density pro- 
files. The King profile provides a good representation of 
the radial distribution of stars in all clusters in the 3Q 
sample, which reinforces their open cluster nature. Since 
it follows from an isothermal (virialized) sphere, the close 
similarity of the radial distribution of stars with a King 
profile suggests that the internal structure of the clusters 
(particularly the core) have already reached some signif- 
icant level of dynamical evolution. We will return to this 
point in Sects. I5"T1 and IB1 The structural parameters are 
given in Table [3 

We quantify the cluster/background star-density con- 
trast by means of the parameter S c (column 4), defined 
as the ratio of the central density (cr(0)) to that of the 
background (ob g ). The S c values in Table [3 can be com- 
pared to those of the rather populous and high-Galactic 
latitude open clusters (Bonatto & Bica l2UU5jl NGC 188 
(b = +22.39°; d - 1.7kpc; M ~ 3800 M ) with S c = 
21.5 ± 2.9, and M67 (b = +31.89°; d Q - 0.9kpc; M ~ 
990 Mq) with S c = 33.7 ± 5.7, and that of the populous 
and rather low-latitude cluster NGC 2477 (b = -5.82°; 



Table 2. 3Q cluster structural parameters 



Object 


o"b g 


&0K 


5 C 


Rcoro 


Rlim 




(stars pc -2 ) 


(stars pc — 2 ) 




(PC) 


(PC) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


Haf9 


4.79 ± 0.08 


32.2 ± 9.2 


7.8 ± 1.9 


0.42 ± 0.10 


4.7 ± 0.6 


Haf 11 


0.54 ± 0.01 


4.0 ± 1.7 


8.6 ± 3.3 


1.09 ± 0.33 


12.1 ± 1.5 


Cz32 


0.55 ± 0.01 


8.4 ± 3.5 


10.4 ± 3.9 


0.81 ± 0.23 


7.8 ± 0.6 


Cz31 


2.75 ± 0.01 


15.1 ± 6.5 


6.4 ± 2.3 


0.53 ± 0.18 


4.0 ± 0.3 


Trl3 


1.92 ± 0.05 


19.0 ± 7.7 


11.2 ± 4.2 


0.59 ± 0.17 


7.1 ± 0.7 



Table Notes. The background stellar density (col. 2) was measured in 
the outermost ring. The central density (col. 3) and core radius 
(col. 5) were derived from the King fit, er(r) — <7"b g + ctok/(1 + 
(r/R core ) 2 ). Col. 4: density contrast parameter 5 C — <t(0)/o , \ ds — 
1 + croK/Vbg- The limiting radius (col. 6) was estimated directly 
from the density profiles. 

d - 1.2 kpc; M - 5300 M s ) with 5 C = 19.7+ 1.9. The 
determination of the density contrast parameter seems to 
depend more on the Galactic latitude than cluster mass. 
Because the 3Q sample is composed of low-latitude clus- 
ters their contrast parameters are significantly lower than 
those of the comparison clusters but higher than the ob- 
served ones (S' c , col. 8 of Table^). As a consequence, some 
underestimation of the limiting radii can be expected in 
most of the clusters studied in the present paper. Deeper 
wide-field photometry would be important to check this 
possibility. 
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Fig. 8. Radial density profiles of the 3Q clusters. The re- 
spective two-parameter King model fits are shown. To al- 
low absolute comparison between clusters, all profiles are 
plotted in the same scale with parsecs in the abscissa and 
stars pc -2 in the ordinate. 



5. Mass functions 

Based on the King profile fits we decided to derive MFs 
(4>{m) = 4^) for the core, halo and overall regions. As 
offset field we consider the region from w 1.5 x Ri im to 
Roxt (Tabled, which provides statistical representativity 
in terms of background distribution of stellar luminosities 
(and mass), as shown in Bonatto, Bica & Santos Jr. ( 2005 ) 
and Bonatto & Bica H2()()5jl . 

In the case of low-latitude clusters it is essential that 
the Galactic field contamination of the CMDs is properly 
taken into account in order to derive the intrinsic lumi- 
nosity and mass distributions of the member stars. To do 
this we first apply the colour-magnitude filter (right pan- 
els of Figs. H to EJ to both cluster and offset field CMDs. 
The filtering process takes into account most of the back- 
ground, leaving a residual contamination. We deal with 
this residual contamination statistically by building the 
luminosity functions (LFs) for each cluster region and off- 
set field. The three 2MASS bands are treated indepen- 
dently, taking into account the 99.9% PSC Completeness 



Limit 6 . Consequently, the faint magnitude limit of each 
LF is J = 15.8, H = 15.1 and K s = 14.3, respectively. We 
take the turnoff as the bright limit to avoid inconsistencies 
in the mass-luminosity relation. For each 2MASS band we 
build a LF by counting stars in magnitude bins from the 
respective faint magnitude limit to the turnoff, both for 
each cluster region and offset field. Considering that the 
solid angle of the offset field is different from that of a given 
cluster region, we multiply the offset field LF by a numer- 
ical factor so that the solid angles match. The intrinsic 
LF of each cluster region is obtained by subtracting the 
respective (i.e. solid angle-corrected) offset-field LF from 
that of the cluster region. Finally, the intrinsic LFs are 
transformed into MFs using the mass-luminosity relation 
obtained from the respective Padova isochrone and dis- 
tance modulus (Sect. EJ- These procedures are repeated 
independently for the three 2MASS bands (Bonatto & 
Bica 2003- The final MF of a given cluster region is pro- 
duced by combining the J, H and Ks MFs into a single 
MF. The resulting core and overall MFs of the 3Q sample 
clusters are shown in Fig. 

We provide in Table|3]parameters derived from the LFs 
and MFs for the core (R < R CO re), halo (R CO ro < R < Riim) 
and overall (core + halo) regions. The number of evolved 
stars (column 3) is calculated by integrating the intrinsic 
LFs for magnitudes brighter than the turnoff. Multiplying 
this by the mass at the turnoff gives an estimate of the 
evolved-star mass (column 4). This procedure produces a 
realistic value of the number of member evolved stars be- 
cause the background contamination was statistically sub- 
tracted from the LF. The observed main sequence (MS) 
mass range is in column 5. The MF slope of the MS stars 
is in column 6. The core MFs are significantly flatter than 
halo MFs (see also Sect. 15. If) . The number of MS stars 
and corresponding mass are derived by integrating the MF 
from the faint magnitude limit to the turnoff. We add to 
these the corresponding values of the number and mass of 
evolved stars to derive the total number of observed stars 
(column 7) and observed mass (column 8). 

The large distances from the Sun of the clusters in 
the 3Q sample preclude detection of sub-Solar mass stars. 
However, an estimate of the total mass locked up in MS 
stars can be made by taking into account all stars from 
the turnoff down to the H-burning mass limit, 0.08 M©, 
assuming the universal IMF of Kroupa l)2()()l|l . in which 
X = 0.3 ± 0.5 for the range 0.08 < m(M©) < 0.5 and 
X = 1.3 ±0.3 for 0.5 < m(M ) < 1.0. For the cases where 
the present value of x is smaller or of the same order than 
Kroupa's we extrapolate our MFs. The resulting extrapo- 
lated values of the number of stars and extrapolated mass 
(added to the corresponding values for the evolved stars) 
are given respectively in cols. 9 and 10 of Tablc|21 We also 
give in cols. 11 and 12 the projected and volume mass 
densities, respectively. 



Corresponding to the Level 1 Requirement, according to 

http:/ /www.ipac. caltech. edu/2mass /releases /allsky /doc/ 'sec6=5al .html 
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Table 3. Measurements and parameters derived from the MFs. 



Evolved MS Observed+Evolved Extrapolated+Evolved 



or 


1 < omirvn 


N* 


mcvol 


Am 


Xms 


N* 


m obs 


N* 


nitot 




P 






(Stars) 


(Mnl 


(Mnl 




(10 2 stars) 


(10 2 M n ) 


(10 2 stars) 


(10 2 M n ) 


(Mn PC~ 2 ) 


(Mn PC~ 3 ) 


1 1 "1 


(91 






C5 1 ) 






\°) 








(121 




Core 


n 

u 


n 
u 


i o— ^ ^ 




u.o in u. i 


o 7 4- n 9 


Q 4 9 


n q 4- n 9 


1 fizt 4 "}0 




Haf9 


Halo 


3 ± 1 


12 ± 5 


1.0-4.1 


±1.35 ± 0.14 


1.6 ± 0.2 


3.0 ± 0.4 


21 ± 15 


8.7 ± 2.9 


13 ± 4 


2.1 ± 0.7 




V../ V Y>1 ctll 


3 ± 1 


12 ± 5 


1.0-4.1 


-1-1 29 ± 1 5 


1.6 ± 0.2 


3.1 ± 0.5 


20 ± 14 


8.5 ± 2.7 


12 ± 4 


2.0 ± 0.6 






6 ± 2 


13 ± 4 


i 7_9 n 




0.4 ±0.1 


\J. i 3Z U. J- 


(17401 

U. t 3Z U. -L 


1.2 ± 0.2 


^9 4 ^ 

t>i 31 (J 


22 4 3 


Haf 11 


Halo 


52 ± 20 


28 ± 9 


1.7-2.1 


-3.25 ± 0.96 


1.3 ± 0.7 


2.6 ± 1.4 


2.1 ± 0.9 


3.6 ± 1.6 


0.8 ± 0.3 


0.05 ± 0.02 




O VPffl 11 
V../ V CI till 


56 ± 20 


28 ± 9 


1.7-2.1 


— 3 61 ±0 82 


1.3 ± 0.7 


2.7 ± 1.3 


2.0 ± 0.8 


3.6 ± 1.4 


0.8 ± 0.3 


05 ± 02 






1 ± 1 


2 ± 2 


i ^—9 n 


_o on _i_ n 99 


n 1 4 n 1 

O.O 3Z U. -L 


n ^ 4- n i 

XJ.O 3Z U. J- 


n ^ + n i 


8 4 1 


OO 3Z O 


OO ZC ij 


Cz 32 


Halo 


11 ±2 


22 ±4 


1.5-2.0 


-1.25 ± 0.74 


0.9 ±0.5 


1.6 ±0.8 


2.9 ±0.9 


3.3 ± 1.1 


1.7 ±0.6 


0.16 ±0.05 




Overall 


11 ±2 


22 ±4 


1.5-2.0 


-1.00 ±0.70 


1.0 ±0.5 


1.7 ±0.8 


3.6 ± 1.2 


3.8 ± 1.1 


2.0 ±0.6 


0.19 ±0.06 




Core 








1.1-2.3 


±0.32 ± 0.40 


0.3 ±0.1 


0.6 ±0.2 


1.8 ±0.8 


1.0 ±0.2 


116 ±28 


165 ± 40 


Cz31 


Halo 


7± 2 


20 ±6 


1.1-2.9 


±1.73 ±0.33 


0.7 ±0.2 


1.2 ±0.3 


11 ± 8.2 


4.1 ± 1.6 


8.2 ±3.1 


1.5 ±0.6 




Overall 


7± 2 


20 ±6 


1.1-2.9 


±1.48 ±0.31 


0.8 ±0.2 


1.3 ±0.3 


11 ±7.9 


4.1 ± 1.5 


8.1 ±3.0 


1.5 ±0.5 




Core 


2 ± 1 


6±3 


1.3-2.8 


-0.91 ±0.37 


0.3 ±0.1 


0.7 ±0.2 


0.5 ±0.1 


0.8 ±0.2 


72 ±23 


91 ±29 


Trl3 


Halo 


23 ±7 


71 ±22 


1.2-3.1 


±1.93 ±0.23 


2.9 ±0.5 


1.4 ±5.3 


58 ±43 


22 ±8.1 


14 ±5 


1.4 ±0.5 




Overall 


24 ±8 


74 ±25 


1.2-3.1 


±1.83 ±0.23 


3.0 ±0.5 


1.2 ± 1.0 


58 ±43 


29 ±8.1 


18 ±5 


1.9 ±0.5 



Table Notes. The number and mass of the evolved stars in cols. 3 and 4, respectively were derived from the LFs. Col. 5: MS 



mass range over which the MF was fitted. Col. 6: MF slope in the MS. 
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Fig. 9. Core and overall mass functions of the 3Q clusters. 
Each panel contains MFs derived from the J, H and Ks 
2MASS photometry. MF fits {<j){m) oc m~( 1+x '>) are shown 
as dashed lines, and the respective MF slopes are given. 
In all cases the core MF is flatter than the overall MF. 



5.1. Considerations on dynamical states 



Mass segregation in a star cluster scales with the relax- 



ation time, defined as t re i a x = 



3 , where t c 



81nN* "cross, 

R/cr v is the crossing time, N* is the total number of 
stars (col. 9 of Table |3J) and er v is the velocity disper- 
sion (Binney & Tremaine I1987fl . t rc i ax is the character- 
istic time scale in which a cluster reaches some level of 
kinetic energy equipartition with massive stars sinking to 
the core and low-mass stars being transferred to the halo. 
We assume a typical velocity dispersion of er v « 3kms _1 
(Binney & Merrifield 1998). Core and overall relaxation 
times are given in cols. 3 and 5 of Table 01 respectively. 
However, since the dynamical evolution depends strongly 
on age and cluster mass (Bonatto & Bica 121)1)51 and refer- 
ences therein) , it is better characterized by the dynamical- 
evolution parameter r = age/t ro i ax . The presence of mass 
segregation and consequently some degree of MF slope 
flattening in the cores of the clusters in the 3Q isamplc 
consistent with the large values of r(core) (Bonatto & Bica 
I2TTO . 

Considering the slope uncertainties, the overall MFs of 
Haf 9 and Cz31 are similar to a standard Salpeter (1955) 
IMF (x = 1.35), while the MF of Tr 13 is steeper. The flat 
MFs of Haf 11 and Cz 32 may be affected by the restricted 
observed MS mass range (AniMs ~ 0.3 M Q and 0.5 M , 
respectively) owing to the large distances from the Sun 
(5.2 kpc and 4.0 kpc, respectively). In all cases the core 
MF resulted significantly flatter than that in the halo, 
which reflects the effects of large-scale mass segregation. 
All these pieces of evidence point to advanced dynamical 
states, especially in the cores of the 3Q sample clusters. 
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200 





Core 


Overall 


Object 




r 


t relax 


T 




(Myr) 




(Myr) 




(1) 


(2) 


(3) 


(4) 


(5) 


Haf9 


0.34 ±0.10 


407 ± 129 


550 ± 31 


2.8 ± 1.2 


Hafll 


0.71 ±0.23 


1250 ± 453 


18 ±6 


49 ± 19 


Cz 32 


0.42 ±0.12 


2660 ± 921 


19 ±5 


58 ±20 


Cz31 


0.73 ±0.37 


248 ± 130 


25 ± 16 


7.3 ±4.8 


Trl3 


0.30 ±0.10 


1060 ± 390 


189 ± 126 


1.7 ± 1.1 



Table Notes. The dynamical-evolution parameter in cols. 3 
and 5 is defined as r = age/t rc i ax . 

6. Comparison with open clusters in different 
dynamical-evolution states 

We compare the 3Q sample clusters with the 631 WEBDA 
open clusters with parameters (Sect, f^) in terms of 
Galactic longitude and latitude, distance from the Sun 
and age. A similar analysis of 13 open clusters studied 
by means of integrated spectra, using reddening and age 
histograms was carried out by Ahumada et al. (2001). 

Galactic longitude and latitude distributions of the 
WEBDA sample are shown in panels (a) and (b) of Fig. 1 101 
respectively. The clusters are rather uniformly distributed 
in I with a slight excess towards the third quadrant, where 
the present clusters are located. This excess probably re- 
flects an observational bias in the sense that low-contrast 
clusters are more easily detected towards regions avoiding 
the Galactic center and bulge. Alternatively, the excess 
may reflect a higher probability of cluster dissolution to- 
wards the Galactic center (e.g. Bergond, Leon, & Guibert 
IMflTl . In Galactic latitude the majority of the WEBDA 
clusters are tightly projected close the Galactic disk, as 
does also the 3Q sample. 

The distribution of the WEBDA open clusters with 
known distance from the Sun is shown in panel (c) of 
Fig. E3 where the distribution peaks at d© « 1.15 kpc. 
Differently from the WEBDA sample, the 3Q sample clus- 
ters distribute rather uniformly in the range d0 ~ 1.9 kpc 
to d Q ~ 5.2 kpc. 

The age distribution of the WEBDA open clusters 
(panel (d)) presents 2 peaks at « 10 Myr and « 100 Myr. 
The 3Q sample clusters distribute rather uniformly in age. 

We conclude that the 3Q sample is not atypical com- 
pared to the WEBDA sample. 

The parameters derived in the previous sections are 
used to check how the 3Q sample clusters fit in the con- 
text of nearby open clusters characterized by different ages 
and dynamical states. Bonatto & Bica ( 2005) presented a 
systematic analysis of a set of 11 nearby open clusters 
with distances from the Sun in the range 0.4 - 1.7 kpc, 
ages from 70 Myr to 7 Gyr and total masses from 400 to 
5 3OOM . These nearby clusters are M26 (NGC6694), 
NGC2516, NGC2287 (M41), M48 (NGC2548), M93 
(NGC2447), NGC5822, NGC2477, NGC3680, IC4651, 
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Fig. 10. Distribution of WEBDA clusters in Galactic lon- 
gitude (panel (a)), latitude (panel (b)), distance from the 
Sun (panel (c)) and cluster age (panel (d)). Filled his- 
tograms: 3Q sample clusters. For clarity the number of 
clusters in panel (a) was multiplied by 4. 



M67 (NGC2682) and NGC 188. The relatively high- 
Galactic latitude and rather populous nature of these 
clusters make them potential sources of accurate intrin- 
sic open cluster parameters. As a result, a set of uniform 
parameters related to the structure (core and overall radii, 
mass and density), dynamical state (core and overall MF 
slopes, dynamical-evolution parameter r = age/t re i ax ), as 
well as age and Galactocentric distance of open clusters 
was obtained. Some correlations among these parameters 
were verified, and a separation of massive (m > 1 000 M Q ) 
and less- massive (m < 1 000 M Q ) clusters was observed in 
some diagrams. The methodology used to derive param- 
eters of the 3Q sample clusters follows the same lines as 
that used in Bonatto & Bica H2005[) . 

In Fig. 1111 we compare the 3Q sample clusters with 
those in the nearby cluster set in terms of structural and 
dynamical-evolution parameters. We include in panels (a), 
(b), (d) and (f) the least-squares fits (and the la borders) 
of parameter relations derived for the nearby clusters in 
Bonatto & Bica lj2UU5|l . 

Within uncertainties the 3Q sample clusters fit well in 
the relation involving core and limiting radii (panel (a)). 
In this case the least-squares fit is given by Rn m = (1.05 ± 
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Fig. 11. Relations involving structural and dynamical- 
evolution parameters of open clusters. Circles: nearby 
open clusters. Triangles: 3Q sample. Filled symbols: clus- 
ters more massive than 1000M Q . Open symbols, m < 
IOOOMq. Dashed lines: least-squares fits to the nearby 
clusters (see text). Shaded areas: la borders of the least- 
squares fits. 



0.45) + (7.73 ± 0.66) x R CO r C , with a correlation coefficient 
CC = 0.95, using the reference sample. 

In the relation involving core and overall mass (panel 

(b) ) the 3Q sample clusters fall slightly off the la border 
of the comparison sample relation, except for Haf 9 which 
fits in it. The least-squares fit is M corc = (14.17 ± 10.02) + 
(0.08 ± 0.01) X Moverall, with CC = 0.92. 

Within the uncertainties the separation observed 
among the massive and less-massive nearby clusters (panel 

(c) ) can be extended to the 3Q sample clusters. 

In the Xcorc x Xovcraii plot (panel (d)) Tr 13, and es- 
pecially Haf 11, deviate from the relation. In this case the 
least-squares fit is Xovcraii = (1-22 ± 0.07) + (0.49 ± 0.08) x 
Xcorc, with CC = 0.90. 

Compared to the nearby clusters of the same age range 
Haf 1 1 seems to present an exceedingly low overall density 
(panel (e)), which might indicate an underestimation of 
the overall mass. This effect may be accounted for obser- 
vationally by the relatively large distance from the Sun 
and low latitude of Haf 11 (Sect. IS~T)) . 

Finally, with respect to the core MF slope vs. 
dynamical-evolution parameter (panel (f)) the least- 
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Fig. 12. Relations of overall density (top panel) and lim- 
iting radius (bottom panel) with the distance to the 
Galactic plane. Symbols as in Fig. 1121 



squares fit is Xcorc = (1.17 ± 0.23) — (3.0 ± 
0.7)exp-( 43 ^ r 1 c 56 ), with CC = 0.82. Compared to the 
nearby cluster set, the cluster cores in the 3Q sample 
present evidence of important dynamical evolution, such 
as flat MF slopes and large values of T corc . This hap- 
pens despite the young age of Haf 9 (~ 140 Myr), Cz31 
(~ 180 Myr) and Trl3 (~ 320 Myr) as compared to the 
nearby open cluster sample. A possible explanation for 
this evidence of advanced dynamical evolution for open 
clusters with ages in the range 100 to 400 Myr may be 
orbital planes closer to the Galactic plane. 

The dependence of limiting radius on distance to the 
plane (|z|) occurs for the less-massive clusters, both for the 
nearby clusters and those in the 3Q sample (bottom panel 
of Fig. lT2*|l . For the massive clusters the dependence of Rii, n 
on \z\ is less important. In the top panel of Fig. lT2l we show 
the relation of the overall density with \z\. A correlation 
shows up for the less-massive clusters in the sense that 
higher-|z| clusters have lower overall densities. However, 
we note that the limiting radius may be underestimated 
(Sect.HJ. The plots in Fig. 1111 fand other parameter rela- 
tions in Bonatto & Bica 2005 ) appear to be useful tools in 
the interpretation of the structure and dynamical states 
of open clusters. 
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7. Concluding remarks 

In this paper we analyzed 5 low-contrast (observed con- 
trast parameter S' c = 1.9 — 4.7) open clusters located 
in the third quadrant by means of 2MASS photome- 
try. The clusters did not have available parameters in 
the literature. The present approach involves Colour- 
Magnitude Diagrams, statistical field-star subtraction, 
colour-magnitude filters, radial density profiles and lumi- 
nosity/mass functions built for the core, halo and overall 
cluster regions. In all cases the CMD morphology result- 
ing from the field-star statistical subtraction is unambigu- 
ously characteristic of open clusters of different ages with 
a well-defined main sequence, turnoff and, for Haf 11 and 
Cz32, the presence of a conspicuous giant clump. 

We derived reddening values E(B — V) in the range 
0.06 to 0.66, ages from 140 Myr to 1.1 Gyr, distances from 
the Sun from 1.9 to 5.2 kpc, and total masses from 360 
to 2 9OOM . 4 out of the 5 clusters are located up to ~ 
3.4 kpc outside the Solar circle. 

By applying the colour-magnitude filter the two- 
parameter King profile could be fitted to the radial den- 
sity distribution of all clusters producing core radii in the 
range 0.42 < R CO ro(pc) < 1.09 and limiting radii in the 
range 4.0 < Rn m (pc) < 12.1. The resulting contrast pa- 
rameters (after applying the colour-magnitude filter) in- 
crease to the range <5 C = 6.4 — 11.2, about 3 times as large 
as the observed ones. 

In all cases the clusters in the 3Q sample present signs 
of advanced dynamical evolution, such as flat core mass 
function slopes, large values of the dynamical-evolution 
parameter, and large-scale mass segregation. This effect 
occurs despite the relatively young cluster ages in part of 
the sample (140 to 320 Myr), resembling dynamical states 
of the nearby Gyr-old clusters (Bonatto & Bica ED05) . 
Proximity to the Galactic plane may have accelerated 
their dynamical evolution. 

We conclude that 2MASS photometry can be used to 
study the large number of unexplored low-contrast open 
clusters when associated with statistical field-star decon- 
tamination (for CMD morphology and isochrone fit) and 
colour-magnitude filters (for structure and mass-function 
analyses) . 
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